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The amplitude of localized quantum states in random or disordered media may exhibit long-range 
exponential decay. We present here a theory that unveils the existence of an effective potential which 
finely governs the confinement of these states. In this picture, the boundaries of the localization 
subregions for low energy eigenfunctions correspond to the barriers of this effective potential, and the 
long-range exponential decay characteristic of Anderson localization is explained as the consequence 
of multiple tunneling in the dense network of barriers created by this effective potential. Finally, we 
show that the Weyl’s formula based on this potential turns out to be a remarkable approximation 
of the density of states for a large variety of one-dimensional systems, periodic or random. 
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Despite more than 50 years of research, many ques¬ 
tions on the exact mechanism of Anderson localization 
still remain open m- One of the most puzzling aspects 
of this phenomenon is the strong spatial confinement of 
the one-particle quantum states, i.e., the exponential de¬ 
cay of the wave amplitude at long range in the absence of 
any confining potential [SHSl- In the interpretation due to 
Anderson, this decay comes from the destructive inter¬ 
ferences between waves traveling from an initial source 
along different propagation pathways in the disordered 
potential um- At distances from the origin much larger 
than the correlation length, the waves statistically almost 
cancel each other out, leading to an exponential decay of 
the amplitude. However, this statistical result says noth¬ 
ing about the detailed nature of the decay. Does it occur 
smoothly, in a continuous way, when the distance goes to 
infinity? Or are there specific places where a transition 
between constructive and destructive interference can be 
observed? 

A recent theory has shown that the precise spatial lo¬ 
cation of such quantum states in a potential V{r) can 
be predicted using the solution u{f) of a simple asso¬ 
ciated Dirichlet problem, called the localization land¬ 
scape HU. While quantum states and their energies are, 
respectively, the eigenfunctions and the eigenvalues of the 

- h? 

Hamiltonian of the system defined as F? =-A + V, 

2m 

the landscape u is defined as the solution of 

Hu = - AuPVu = l , (1) 

2m 

the boundary conditions being either Dirichlet, Neu¬ 
mann, or periodic. In this theory, the localization sub- 
regions are delimited by the valley lines of the graph 
of u. This property directly derives from a fundamen¬ 
tal inequality satisfied by any eigenfunction oi H with 
eigenvalue E, normalized so that its maximum amplitude 


is equal to I: 


< E u{r) . (2) 


In other words, the small values of u along its valley lines 
constrain the amplitude of to be small along the same 
lines and, as a consequence, localize low energy eigen¬ 
functions inside the regions enclosed by these lines m- 
We unveil here a different—and much more powerful 
role—played by u, by showing that the function W = Iju 
can in fact be interpreted as a confining potential that 
is responsible for the exponential decay of the Anderson 
localized states even far from its main localization sub- 
region. To that end, the original Schrodinger equation is 
transformed by introducing an auxiliary function (p such 
that = uip. Expressing that is an eigenvector of the 
Hamiltonian leads to 

^ + {uip) = E Uip . (3) 

Developing this equation and accounting for the defini¬ 
tion of u in Eq. 0 gives 

h? A 1 rn /.N 

— -— Ap) — 2- -• V(/9 H— (^ = E(p . (4) 

2m 2m u u 

The additional first order term proportional to V(^ can be 
inserted into the second order term which finally yields: 
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^ div 


EWif 


Eip . 


( 5 ) 


One can see that the auxiliary function cp = fj/u thus 
obeys a Schrodinger-type equation in which the original 
potential V (r) has disappeared. Instead, a new function 
W(r) now plays the role of “effective confining potential.” 
One first notices that, since u is a solution of Eq. 0, 
W is indeed homogeneous to an energy. Moreover, the 
valleys of u which are the boundaries of the localization 
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FIG. 1. (a) Random piecewise constant potential, (b) Effec¬ 
tive confining potential W computed by solving Eq. then 
taking the reciprocal of u. (c) Probability amplitude of the 
fourth state j'lps) of energy E 3 = 1.88 in a logarithmic scale 
superimposed on the effective potential W (the dashed line). 
The wells of W for this state are defined by the intersec¬ 
tions with the horizontal line W — and are outlined in 
grey. One can see that the decay of the eigenfunction oc¬ 
curs mostly outside these wells, i.e., in the barriers of VF, 
and stops inside these wells, (d) Amplitudes of three differ¬ 
ent quantum states (V^o, and in a logarithmic scale, 
superimposed with the estimates obtained using the Agmon 
distance (the dashed lines). These estimates are equal to the 
right-hand side of Eq. The Agmon estimates follow in de¬ 
tail the decay of the actual amplitudes down to values smaller 
than 10 “^. 

subregions m also correspond to the crest lines of this 
new potential. These crest lines act as barriers for the 
auxiliary function ip. 

The demonstration that W plays the role of an effective 
potential derives from the following equality satisfied by 
any quantum state l^p): 

= {upp (E^) + (V’l^it/’). (6) 

The detailed proof of this equality is given in the Supple¬ 
mental Material m- It not only shows that, regardless 
of its kinetic energy, the energy of a quantum state \ 2 p) 
can never be smaller than the one it would have in a po¬ 
tential W (r), but also shows that the difference (W — E) 
can be used to build an Agmon distance pE (Fi) that con¬ 
trols the decay of %jj{r) in the regions where < FF, as 
dictated by Agmon’s inequality [nun]. This distance is 
defined as 

pE (n, r5) = min (^ J {W {f) - E)^ (7) 

where the minimum is computed on all paths 7 going 
from rl to r 2 . The control on the amplitude '0(r) of 


an eigenfunction centered in ro of energy E is expressed 
through the inequality 

^ . ( 8 ) 

This exponential decay can be observed in Fig. which 
displays the localization of quantum states in a random 
potential with periodic boundary conditions. The poten¬ 
tial is piecewise constant on intervals of length 1, with 
values following a uniform law between 0 and FFnax = 8, 
where units of Y? j2m are considered [see Fig. Ba)]. The 
domain length E = 64 thus corresponds to the total 
number of such intervals. Figure [^b) displays the effec¬ 
tive confining potential FF = 1/rt computed from solv¬ 
ing Eq. Q. This effective potential (the dashed line) is 
superimposed in Fig.j^c) with the logarithmic plot of the 
third excited state (V^s). In a log scale, an exponential 
decay of the state amplitude translates into a linear drop. 
The wells of FF for are the locations where FF < E3, 
and they are outlined in grey. One can observe that the 
decay of -03 occurs exactly at the barriers of FF and stops 
across the wells. Fig. Bd) presents a comparison between 
the amplitudes of 7/^0 5 ^.nd 7/^3 (the solid lines) and 
the estimates obtained using the Agmon distance build 
from FF (dashed lines). We find that the amplitudes and 
their estimates are almost identical, and this result is ro¬ 
bust for higher energies and for many realizations of the 
random potential. 

The confining properties oi 1/u are even more inter¬ 
esting in two dimensions. In the case of a 2D random 
potential where no clear localization region can be out¬ 
lined, it has already been observed that the localization 
landscape u exhibits marked valleys that determine the 
localization subregions m- Here we show that 1/u expo¬ 
nentially controls the confinement of the quantum states 
in the entire domain through the heights and widths of 
its barriers. This is of particular interest in cases where 
the semiclassical approach fails, as, for instance, for a 
Boolean-type potential (a random potential that can take 
only two values). Figure [^a) displays a realization of a 
2D Boolean potential, i.e., a potential that can take only 
two values, FFnin = 0 or FFnax = 4. In the present com¬ 
putation, the domain is divided into 40 x 40 small unit 
squares, and the potential is piecewise constant on each of 
these unit squares, taking the value FFnin with probabil¬ 
ity Po = 0-6 and the value FFnax with probability p = 0.4 
so that the V = 0 region has a large chance to percolate 
throughout the domain, as happens for the specific real¬ 
ization. Consequently, a classical particle would never be 
confined by such a potential, independent of its energy. 

Both the effective potential FF and the fundamental 
quantum state are obtained using 320,000 triangular 
Lagrange finite elements of degree 3, with a code based 
on the FEniCS finite element software environment m- 
Figure [^b) displays a color representation of FF, and 
its crest lines are computed using a watershed algorithm. 
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(c) (d) 


FIG. 2 . (a) Boolean potential V. The black region corre¬ 

sponds to y = Vmax = 4. The white region, corresponding to 
y = 0 , occupies about 60% of the domain and percolates from 
the center to the outer boundary of the domain, (b) Color 
representation of VF = 1/u, the effective potential as defined 
in Eq. ([5), superimposed by the crest lines of this potential 
computed using a watershed algorithm, (c) Color representa¬ 
tion of log 3 ^Q[' 0 o(T)], 'tpo being the fundamental eigenfunction 
of energy Eq = 0.503, superimposed with the crest lines of W. 
(d) The logarithmic gradient of 'ipo, thresholded so that values 
> 1 are in black, and superimposed with the crest lines of W. 
The clear match indicates that the exponential decay of '0 oc¬ 
curs exactly at the boundaries of the localization subregions 
determined by the effective potential W. 


Figure]^ c) presents the amplitude of overlapped with 
the same lines. One can notice how the decay of the am¬ 
plitude closely follows the subregions delimited by these 
lines. 

In order to examine in detail the long-range decay of 
this quantum state, the quantity |Vlog(|'0|)| which can 
be interpreted as the inverse of the local length of decay 
of the wave function, is plotted in Figure |^d). Superim¬ 
posing the level sets of this quantity with the crest lines 
of 1/u (which are also the valleys of the landscape u) re¬ 
veals a very strong correlation: the local length of decay 
increases substantially in the vicinity of the barriers of W. 
The long-range exponential decay in Anderson localiza¬ 
tion does not occur uniformly but rather appears as the 
consequence of successive and cumulative decays across 
the dense network of barriers generated by W. This be¬ 
havior can be observed consistently in many trials and 
for all localized eigenfunctions. 

In summary, the decay of the original eigenfunction 
ip = u (p away from its maximum thus originates from 
two concurring contributions: 

(a) First, u is small near the valleys surrounding the 


maximum of -0, reducing accordingly the amplitude 
of ip wherever E u{r) <1 (see Ref. [H]). 

(b) Second, p decays through the barriers of FF, which 
can be interpreted as “quantum tunneling” of the 
auxiliary function when FF is larger than E. 


So, not only does the amplitude of the quantum state pj 
decrease in the vicinity of the valleys of the landscape 
but it also decays when crossing the faraway valleys, as 
long as the value of FF is larger than the energy E. This 
decay is directly governed by the exponential of the Ag- 
mon distance build from FF. 

One can therefore conclude that the complicated in¬ 
terferential pattern of wave localization in V is trans¬ 
lated into the picture of a classical confinement, observed 
through the glass of the wells and the barriers of the new 
effective potential FF. It is thus possible to graphically 
identify the localization subregions at energy E just by 
Hooding the effective potential FF up to the height E^ 
and by then observing the extent of the flooded basins. 
Each basin can be considered a local oscillator, the entire 
system then appearing as a set of oscillators, each inde¬ 
pendent at low energy and weakly coupled by quantum 
tunneling through the boundaries of the basins. One has 
to underline that for energies higher than the maximum 
of FF, there is no more classical confinement by FF, and 
the state localization in ID or 2D comes from the ran¬ 
domness present in the effective potential. 

It is interesting to note that FF can also be interpreted 
as a regularized version of the original potential V. This 
can be seen by rewriting Eq. as 


V-W = 


2m u 


(9) 


However, unlike in a classical smoothing procedure, the 
smoothing scale here is not constant but varies spatially 
depending on the value of and hence of V. The 
smoothed effective potential FF resulting from this non¬ 
linear operation favors the emergence of well-formed wells 
surrounded by barriers, even in situations where none are 
visible in the original potential. In that sense, the behav¬ 
ior of FF is much closer to that of a classical confining 
potential than V. 

The emergence of localized states triggered by the 
quenched disorder may also strongly perturb the one- 
electron density of states [TTHIS]. In some cases, it may 
shift the energy of the fundamental state, leading to an 
enhancement of the already existing gap in crystalline 
semiconductors m or amorphous semiconductors m or 
the creation of a pseudogap in superconducting materi¬ 
als. Very generally, Weyl’s law states that the integrated 
density of states (IDOS) N{E) at a given energy E can 
be approximated by the volume (properly normalized) in 

the phase space (x^kj that can be explored by a classical 
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particle of mechanical energy smaller than E: 


N{E) 




d^x d^k 


H(x,k)<E 


B^+V(x)<E 


d^x d^k . 



where n is the spatial dimension. For example in the 
local band structure theory of semiconductors, the local 
density of states is obtained by assigning to V the value 
Ec{x)^ the bottom edge of the conduction band. How¬ 
ever, Weyl’s formula is only valid in the asymptotic limit 
E +OC and can be very inaccurate at low energies. 

It follows immediately from Eq. ( [Tq| ) that Weyl’s for¬ 
mula in one dimension writes N{E) ^ Ny^E) where 


(a) (b) 



= IIL 


!q^+V(x)<E 

-il\L 

y/2m 


dx dk 


'^<E-V{x) 


dk 


dx 


7vh 


[ [E-V{x)]^^^dx . 

Jv{x)<E 


( 11 ) 



(d) 
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Because of the analogy with classical mechanics that 
is implicitly behind Weyl’s law, one can easily under¬ 
stand why this approximation is poor in quantum sys¬ 
tems where wave interference plays a major role. This 
can be seen in Fig. where the density of states for three 
typical cases of one-dimensional potential are examined: 
(a) random with uniform law, (c) random boolean, and 
(e) periodic. For each type of potential represented on 
the left, both the true counting function N{E) (the step¬ 
wise solid line) and the approximation Nv{E) (the dot¬ 
ted red curve, obtained from Weyl’s approximation) are 
represented. One can see that in all three cases, these 
two curves differ significantly as expected. 

We now compute a new approximation of the count¬ 
ing function, this time based on the potential W deduced 
from u. The reason for doing so lies in the fact that W 
has been shown to behave as an effective potential energy 
[see Eq. m with identifiable wells and barriers, in the 
spirit of a classical mechanical system. In all three of the 
frames (b), (d), and (f), the smooth curve (the dotted 
red line) is the approximation Nw{E) obtained by in¬ 
serting the effective potential W into Eq. ( [IT] ). One can 
immediately notice the remarkable agreement between 
the original counting function and Nw{E)^ whereas the 
standard approximation Ny using the original potential 
totally fails to predict the density of states. In particular, 
Nw {E) detects precisely the shift of the lower edge of the 
conduction band induced either by disorder [see Fig.j^b)] 
or by the periodicity of the potential [see Fig.j^f)]. This 
lower bound to all energy values corresponds to an added 
gap in the IDOS that is reminiscent of the missing density 
of states at lower energy in Anderson localization [22] . At 


(e) (f) 


FIG. 3. (Left panels) Three different types of potentials: 
(a) random with uniform law on [0 1], (c) random Boolean (0 
or 1), and (e) periodic with 256 periods, F = 0 on the first 
half of the period, F = 1 on the second half. (Right panels) 
For each potential, the counting function N (the solid black 
line) is represented, together with Ny (the blue dash-dotted 
line), Weyl’s approximation using the original potential F [as 
defined in Eq. 0], and Nw (the red dotted line), the same 
with W. Notice in all cases the remarkable agreement be¬ 
tween N and Nw- 


higher energy, all three curves catch up, as they asymp¬ 
totically follow Weyl’s law without potential, which is 
proportional to 

In summary, a new conceptual tool has been intro¬ 
duced for understanding the localization properties of 
quantum states. This object is the reciprocal of the local¬ 
ization landscape introduced in Ref. m and can be inter¬ 
preted as an effective potential. We have shown that this 
effective potential not only determines the boundaries of 
the localization regions but also controls the long-range 
decay of the quantum states through an Agmon metric. 
In particular, the exponential decay observed in Ander¬ 
son localization has been shown not to occur uniformly, 
but rather to be concentrated across the barriers of this 
effective potential. Therefore, the transport between ad¬ 
jacent basins of the effective potential involves an “ef¬ 
fective quantum tunneling.” In short, one can say that 
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the effective potential captures the interference pattern 
created by waves in the original potential and converts 
it into a semiclassical picture of confining potential. Fi¬ 
nally, this property enables us to build an approximation 
of the integrated density of states based on Weyl’s law. 
Applied to various types of ID potentials, either random 
or deterministic, this approximation showed a remarkable 
agreement with the actual IDOS, far more accurate than 
what can be obtained using Weyl’s law, with or without 
the original potential. 

Because it captures the deep localization properties 
of a complex or random potential, the effective confin¬ 
ing potential is a very promising tool for understanding 
the features of quantum waves in disordered media. In¬ 
deed, our procedure reveals the system as a partition of 
weakly coupled oscillators, and it thus efficiently realizes 
an approximate diagonalization of the Hamiltonian for 
the portion of the spectrum spanned by localized eigen¬ 
functions. With this information at hand, we can calcu¬ 
late the carrier distribution in random alloys such as, for 
instance, the active layers of GaN-based Quantum-well 
light-emitting diodes [23]. Finally, our approach brings 
a new perspective to such open questions as the nature 
of the transition from localized to delocalized states in 
three dimensions (the mobility edge) and the onset of 
localization in many body systems. 
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